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We consider spherically symmetric Einstein-massless-scalar field equations with negative cosmo¬ 
logical constant in five dimensions and analyze evolution of small perturbations of anti-de Sitter 
spacetime using the recently proposed resonant approximation. We show that for typical initial 
data the solution of the resonant system develops an oscillatory singularity in finite time. This re¬ 
sult hints at a possible route to establishing instability of AdS under arbitrarily small perturbations. 


Introduction. A few years ago two of us gave numeri¬ 
cal evidence that anti-de Sitter (AdS) spacetime in four 
dimensions is unstable against black hole formation for 
a large class of arbitrarily small perturbations pQ. More 
precisely, we showed that for a perturbation with ampli¬ 
tude e a black hole forms on the timescale 0(£~ 2 ). Using 
nonlinear perturbation analysis we conjectured that the 
instability is due to the turbulent cascade of energy from 
low to high frequencies. This conjecture was extended to 
higher dimensions in [2], 

Since the computational cost of numerical simulations 
rapidly increases with decreasing e, our conjecture was 
based on extrapolation of the observed scaling behavior 
of solutions for small (but not excessively so) amplitudes, 
which left some room for doubts whether the instabil¬ 
ity will persist to arbitrarily small values of e (see e.g. 
0)- To resolve these doubts, in this paper we validate 
and reinforce the above extrapolation with the help of 
a recently proposed resonant approximation 0HB]. The 
key feature of this approximation is that the underlying 
infinite dynamical system (hereafter referred to as the 
resonant system) is scale invariant: if its solution with 
amplitude 1 does something at time t , then the corre¬ 
sponding solution with amplitude e does the same thing 
at time t/e 2 . Moreover, the latter solution remains close 
to the true solution (starting with the same initial data) 
for times < e~ 2 (provided that the errors due to omission 
of higher order terms do not pile up too rapidly). Thus, 
by solving the resonant system we can probe the regime 
of arbitrarily small perturbations (whose outcome of evo¬ 
lution is beyond the possibility of numerical verification). 

For concreteness, in this paper we focus our attention 
on AdSs (the most interesting case from the viewpoint of 
AdS/CFT correspondence); an extension to other dimen¬ 
sions is straightforward and will be presented elsewhere. 

Model. For the reader’s convenience, let us recall from 
ED 0 the general framework for studying the spheri¬ 
cally symmetric scalar perturbations of the AdS space- 
time. The five dimensional asymptotically AdS space- 
times are parametrized by the coordinates (i, x, to) € 
(—00,00) x [0,7t/2) x S 3 and the metric 

ds 2 = ^ (—Ae~ 2S dt 2 + A~ 1 dx 2 + sin 2 xdto 2 ) , (1) 

cos-a; v ' 


where I 2 = —6/A, dui 2 is the round metric on § 3 , and A, 
5 are functions of (t,x). For this ansatz the evolution of 
a self-gravitating massless scalar field </>(t, x ) is governed 
by the following system (using units in which 8irG = 3) 

<9 t $ = d x (Ae~ s n) , c> t n = — \—d x (tan 3 a: Ae~ s <&) , 
v taiU:r 

(2) 

d x A = —7—■———— (1 — A) — sin a; cos a; A ($ 2 + II 2 ) , 
sin x cos x 

(3) 

d x 8 = — sin a; cos a; ($ 2 + II 2 ) , (4) 

where = d x <j> and II = A~ 1 e s dt4>. To ensure smooth¬ 
ness at spatial infinity and finiteness of the total mass M 
we impose the boundary conditions (using p = 7t/2 — x) 

<t>(t,x) = f ao (t)p i + 0(p e ) , S(t,x) = 600(f) + O (p s ) , 
A(t,x) = l-Mp 4 + 0(p 6 ) , (5) 

where the power series expansions are uniquely deter¬ 
mined by M and the functions foo(t), 600(f)- We use the 
normalization 6(t, 0) = 0, hence t is the proper time at 
the center. We will solve this system for small smooth 
perturbations of the AdS solution </> = 0,A = l,<5 = 0. 

Resonant approximation. As follows from equation 
linearized perturbations of AdSs are governed by the op¬ 
erator L = — tan _3 a; 3 x (tan 3 :r d x ). This operator is es¬ 
sentially self-adjoint with respect to the inner product 
(f,g) '■= Jo / \f(x)g(x) tsm 3 xdx. The eigenvalues and 
orthonormal eigenfunctions of L are w 2 = (2 n + 4) 2 
(n = 0,1 ,...) and 

e n (x) = 2 J ^ 1 + + 3 ) cos 4 xP/ 1 ’ 2) (cos(2a:)) , (6) 

V n+1 

where P^’ 2 \x) is the Jacobi polynomial of order n. 

After these preliminaries we are prepared to introduce 
the resonant approximation. To avoid technicalities, let 
us first illustrate this approach in the case of the cubic 
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wave equation on the fixed AdSs background 1 

dttfi + L(f> + sec 2 x (f> 3 = 0. (7) 

Inserting the mode expansion <fi(t,x) = c n(t)e n (x) 
into ([7]) we get an infinite system of coupled oscillators 


d 2 « 


dt 2 


n 2 T 

' ^rS'n — / J Ijkin CjCkQ : 

jkl 


where the coefficients Ijkin = — (e^e/ce; sec 2 x, e n ) deter¬ 
mine interactions between the modes. To factor out fast 
linear oscillations in ([ 8 ]), we change variables using the 
variation of constants (a.k.a.“interaction picture”) 


c n = /3 n e lulr,t + P n e 

^ = *w n (/3 n e 


-iUlnt 

(9) 

- P n e~ iUrit ) • 

( 10 ) 

into 


CjCkCi e ~ lul,lt , 

( 11 ) 


jkl 

where each Cj in the sum is given by (9), thus each term in 
the sum has a factor e _lS4t , where D = oj n ± ojj ± u>k ± uii. 
The terms with Q = 0 correspond to resonant interac¬ 
tions, while those with Q 7 ^ 0 are non-resonant. 

Passing to slow time r = e 2 t and rescaling f) n (t ) = 
ea n (r), we see that for e going to zero the non-resonant 
terms oc e _ * nr / e are highly oscillatory and therefore neg¬ 
ligible (at least for some time). Keeping only the resonant 
terms in 0 (which is equivalent to time-averaging), we 
obtain the infinite autonomous dynamical system (which 
we shall call the resonant system ) 


2 ica, 


dot n 

dr 


^ ( Ijkin OljC^kOll 

jkl 


( 12 ) 


where the summation runs over the set of indices {j, k, 1} 
for which D = 0 and Ijkin 7 ^ 0 (due to the fully resonant 
non-dispersive spectrum of L and vanishing of some coef¬ 
ficients Ijkin, this set reduces to {jkl \j + k — l = n}; see 
footnote 3 in 0)- Note that the system ( ]12j ) is invariant 
under the scaling a n (r) —► e _ 1 a n (r/e 2 ). It is routine 
to show that the solutions of 0 starting from small 
initial data of size e are well approximated by the solu¬ 
tions of (12) on the timescale 0(e~ 2 ) 2 . In other words, 


on this timescale the dynamics of solutions of the cubic 
wave equation ([7j is dominated by resonant interactions. 

For the system (2-4) the derivation of the resonant sys¬ 
tem is similar but technically more intricate because, to 
begin with, one has to integrate out the constraints which 
(at the lowest order) results in nonlocal cubic nonlinear¬ 
ities in the derivatives of <p. Despite these complications, 


( 8 ) the resonant system has the same form as ( 12 ), namely 


2-A’n ^ ' Cjkln OLjOLkOil . 

j-\-k—l—n 


(13) 


except that now the interaction coefficients Cjkin are 
given by much more complicated expressions involving 
integrals of products of eigenfunctions © and their 


derivatives. The system (13) was first derived in [4] and 
[ 1 ] using the multiscale perturbation methods 3 and soon 
afterwards in [ 6 ] using the averaging method. It is sim¬ 
pler than the full system, yet still too difficult to be han¬ 
dled by purely analytic means, hence in what follows we 
analyze it using numerical and asymptotic methods. 

Results. We solved in parallel the full Einstein equa¬ 


tions (2-4) and the resonant system (13) for a variety 
of the same small initial data. To illustrate the results 
(which we believe are universal), we present them for the 
two-mode initial data with energy (almost) equally dis¬ 
tributed among the modes 


0 ( 0 , x) = e ( ^eo(x) + ^ei(a;) 


n(o,x) = o. (14) 


Let us point out that, since small one-mode data and 
their perturbations enjoy quasiperiodic evolutions ms 
0 , the two-mode data are in a sense minimal in what it 
takes to initiate the turbulent cascade (cf.lIUl). 

The numerical simulations of the full Einstein equa¬ 
tions (2-4) were previously reported in [2] in the case 
of Gaussian initial data. The evolution of the two-mode 


data (14) looks similar. For all considered small values of 
e we observe collapse to a black hole in time tu{s) ~ £~ 2 
(see Fig. 1). This scaling suggests that the instability 
should be seen in the resonant approximation. In what 
follows, we confirm this expectation and thereby give sup¬ 
port to the conjecture that the instability is present for 
arbitrarily small perturbations. 

For the numerical computation, the resonant sys¬ 


tem (13) must be truncated at some (possibly large) in¬ 
dex N. As a compromise between the accuracy 4 and the 


1 The study of this equation has been argued (see [7] and references 
therein) to provide insight into the problem of stability of AdS. 
We do not share this view (as explained in 0) and use equation 

o 0 solely for illustrative purposes, not as a toy model. 

2 More precisely, if (3 n (t ) and a n (r) are solutions of ( |11| ) and 
l |12| ), respectively, and (3 n (0 ) = £Q!n(0) for each n, then \/3 n (t) — 
ea n (e 2 t) \ < e for t < e~ 2 . 


3 In the multiscale approach, the system ( |13| ) follows from elimi¬ 
nation of secular terms due to resonances at the third order of 
perturbation expansion. The fact that all secular terms can be 
removed in this way has been sometimes misunderstood as evi¬ 
dence for stability of AdS. 

4 See nn for the analysis of the combined error due to truncation 
and averaging for semilinear wave equations. 
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FIG. 1: Time of horizon formation vs amplitude in the evo¬ 
lution of initial data (14). The solid line depicts the fit of 
the function —2 lne + a + be 2 to the numerical data In trr(e). 
From this fit we obtain th '■= lim £ _>o e= e a ~ 0.514. 


computational cost, we choose here N = 172. To solve 
the truncated resonant system (TRS) numerically we use 
the 6th-order Gauss-Runge-Kutta method 5 . The initial 
data for the TRS corresponding to |l4| ) are 

a 0 (0) = l/8, ai(0) = l/12, (a„(0)) n > 2 = 0. (15) 


To describe and analyze the behavior of solutions, it 
is convenient to use the amplitude-phase representation 
a n = A n e lBrl in terms of which the resonant system (13) 
takes the following form [B] 

= F Sj k i n AjA k Ai sin(B n + B t - B 0 - B k ) 

j -\-k—l—n 
j^n,k^n 


(16) 


2 


dB n 

dr 


— TnA^ + ^ ' RjnA 


j^n 


T A n ^ ) Sjki n AjA k Ai cos ( B n T Bi Bj B k ). 

(17) 


j-\-k—l=n 
j^n,k^n 


where T n — C nnnn , Rj n — C 'njjn ^ j 7 ^ ^5 B-nd Sj^i n — 
Cjkin if b 11 four indices are different. Explicit expressions 
for these coefficients are given in the appendix A in [5]. 

Evolving the data (15), we find that higher modes are 
quickly excited (see Fig. 2). For early times the ampli¬ 
tudes grow at the polynomial rate A n (r) ~ r ra_1 while 
the phases evolve approximately linearly. At a later time 


FIG. 2: The amplitude spectra (for several times r). 
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FIG. 3: Evolution of a sample high mode. 


the frequencies of oscillations begin to grow rapidly (see 
Fig. 3). This highly oscillatory behavior accumulates at 
a finite time causing numerical difficulties. We find that 
the time-step of numerical integration, for which the al¬ 
gorithm is convergent, tends to zero as the cutoff N in¬ 
creases. This suggests that the solution of the resonant 
system (TV = oo) develops an oscillatory singularity in 
some finite time r* 6 . 

To give better evidence for blowup and analyze its 
character we proceed in the spirit of the analyticity strip 
method mm, namely we make the following asymp¬ 
totic ansatz for the amplitudes 

A n (r) ~ for n > 1. (18) 

Fitting this formula to the numerical data we obtain the 
time dependence of the exponent j(t) and the ‘analytic¬ 
ity radius’ p(r). As shown in Fig. 4, it appears that p(r) 


5 Due to the oscillatory character of solutions, explicit schemes 

turn out to be unstable. More details of the numerical method 
and its validation are given in the supplementary material ac¬ 
companying this note. 


For any finite N the solution of TRS cannot blow up and, with 
sufficient resolution, can be numerically continued past r*, how¬ 
ever this ‘afterlife’ is an artifact of truncation. 
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tends to zero in a finite time r* (with po = —p'(r*) > 0 ) 


confirming that the solution of the resonant system (13) 


becomes singular at r*. The fit also reveals that the 
asymptotic power-law amplitude spectrum has the expo¬ 
nent lim-r-^ 7 (r) = 2 . 
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FIG. 4: The radius of analyticity p(r) obtained by fitting the 
formula ( |18| ) to the amplitude spectrum. The point of this 
(notoriously poor) fit is to show that p(r) hits zero in some 
finite time t*, not to determine r* precisely (cf. footnote 6). 


Guided by these numerical findings we will now con¬ 
struct an asymptotic solution of the resonant system that 
becomes singular in finite time. We assume that for large 
n and r —> t* 

A„{t) ~rr 2 e- p ° (T *- T)n . (19) 

To solve for the phases, we note the following asymptotic 
behavior of the interaction coefficients 


T n 


Bjn 


n 2 j\ 




A 4 S. 


jkln • 


( 20 ) 


Notice that the latter implies that 


^ S jkln (jkl)~ 2 = 0( 1), 


( 21 ) 


j-\-k—l=n 

j^n,k^n 


that is the sum does not dependent of n. Plugging (19) 
into Eq.(17) and using (|20| and (21), we see that for 


r —► r* the r.h.s. of Eq.(17) is dominated by the term 


RjnAj 


b „ n 2 ln ( Tst _ r ) . 


( 22 ) 

thus the derivatives blow up logarithmically. More¬ 
over, it follows from the above that for large n the phases 
B n behave linearly with n, hence B n + Bi — Bj — B^ ~ 0 
for the resonant quartets. This implies that both sides 
of Eq.(16) are (approximately) independent of n, reas¬ 
suring that the ansatz (19) is self-consistent. Numerical 


simulations indicate that the asymptotics of blowup just 
described is in fact universal; this is illustrated in Fig. 5 


in the case of two-mode data (15). 



FIG. 5: Evidence for the logarithmic blowup. The solid line 
represents the fit of the theoretical prediction a n ln(r* —r)+b n 
to the numerical data for n = 96. Performing this fit for 
all n > 20 we find, in accord with the asymptotic analysis, 
that the coefficients a n and b n vary linearly with n, while 
r* ss 0.509 does not depend on n. Note that t* agrees very 
well with the time of collapse th given in the caption of Fig. 1. 


Conclusion. To summarize, we have constructed the 
asymptotic solution of the resonant system that becomes 
singular in finite time and gave numerical evidence that 
this solution acts as a universal attractor for blowup. The 
key question is how to transfer this blowup result from 
the resonant system to the full system. On one hand, we 
see that the resonant approximation reproduces the am¬ 
plitudes of true solutions remarkably well almost all the 
way to collapse. This is illustrated in Figs. 6 and 7 where 
we compare the energy spectrum and the growth of the 
Ricci scalar at the origin computed in parallel using the 
full Einstein equations and the resonant system. On the 
other hand, the resonant approximation does not work 
so well for the phases 7 . For this reason (and because of 
possible breakdown of the cubic approximation near col¬ 
lapse), it is not clear to us what (if any) is the physical 
interpretation of the oscillatory singularity for the reso¬ 
nant system 8 . Nonetheless, the fact that solutions of the 
resonant system blow up in finite time (for typical initial 
data) strongly indicates that the corresponding solutions 
of the full system collapse on the timescale 0(e~ 2 ). 
Acknowledgments: This work was supported in part by 
the Polish National Science Centre Grant no. DEC- 
2012/06/A/ST2/00397. The computations were carried 
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7 Note that the highly oscillatory behavior of solutions is in tension 
with the idea of time averaging. 

8 Since the resonant system has no scale, it is tempting to speculate 
that there is a relationship between the oscillatory singularity 
and Choptuik’s discretely self-similar solution m- 
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FIG. 6: Energy spectra for the full system (dotted lines) and 
TRS (solid lines) at sample late times (cf. Fig. 2). The power- 
law spectrum n ~ 2 unfolds as the solutions collapse/blow up. 
Here e s=s 0.0079 (the smallest amplitude used in simulations). 



[1] P. Bizon, A. Rostworowski, Phys. Rev. Lett. 107 , 031102 
(2011), [arXiv: 1104.3702] 

[2] J. Jalmuzna, A. Rostworowski, P. Bizon, Phys. Rev. D 
84, 085021 (2011), [arXiv: 1108.4539] 

[3] F.V. Dimitrakopoulos, B. Freivogel, M. Lippert, I.-S. 
Yang, [arXiv:1410.1880] 

[4] V. Balasubramanian, A. Buchel, S.R. Green, L. Lehner, 
S.L. Liebling, Phys. Rev. Lett. 113, 071601 (2014), 
[arXiv:1403.6471] 

[5] B. Craps, O. Evnin, J. Vanhoof, JHEP 10 , 48 (2014), 
[arXiv:1407.6273] 

[6] B. Craps, O. Evnin, J. Vanhoof, JHEP 01 , 108 (2015), 
[arXiv:1412.3249] 

[7] P. Basu, C. Krishnan, P.N. Bala Subramanian, 
[arXiv:1501.07499] 

[8] http: //www. etc . cam. ac .uk/activities/adsgrav2014/ 
Slides/Slides_Bizon.pdf 

[9] M. Maliborski, A. Rostworowski, Phys. Rev. Lett. Ill, 
051102 (2013), [arXiv: 1403.5434] 

[10] P. Bizon, A. Rostworowski, [arXiv: 1410.2631] 

[11] M.S. Krol, Math. Methods Appl. Sciences 11, 649 (1989) 

[12] C. Sulem, P.-L. Sulem, H. Frisch, J. Comput. Phys. 50, 
138 (1983) 

[13] P. Bizon, J. Jalmuzna, Phys. Rev. Lett. Ill, 041102 
(2013), [arXiv: 1306.0317] 

[14] M.W. Choptuik, Phys. Rev. Lett. 70, 9 (1993) 


FIG. 7: Upper envelope of n 2 (t, 0) in the evolution of initial 
data ( |14| | with e = (27r)~ 3 / 2 2 _p for p = 1,2,3 (dotted lines). 
As £ decreases, the rescaled quantities e~ 2 H 2 (e 2 t, 0) approach 
a limiting curve. The corresponding solutions of TRS (solid 
lines) appear to approach the same curve as N increases. 












